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ABSTRACT 


Various  aspects  of  guidance  aoid  control  for  ballistic 
missiles  are  presented.  A  discussion  of  the  required 
velocity  vector  is  given  with  severed  examples  based 
upon  the  flat  earth  model.  Error  analysis  procedures 
involving  inertial  reference  and  measurement  sys¬ 
tems  are  treated  in  some  detail.  . 

Missile  steering  techniques  are  approached  in  such  a 
way  that  stability  is  not  a  problem.  Procedures  for 
insuring  that  the  missile  velocity  will  be  controlled  to 
the  required  velocity  without  payload  loss  penalties 
are  given.  This  is  accomplished  by  introducing  tra¬ 
jectory-shaping  functions  that  cause  the  guided  bal¬ 
listic  missile  to  closely  follow  the  nominal  reference 
trajectory. 

The  steering  discussion  initially  assumes  apointmass 
missile.  Control  system  -  conipli cations  that  result 
from  actual  missile  dynamics  are  treated  separately. 
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I.  GENERAL  DISCUSSION 


The  purpose  of  a  ballistic  missile  guidance  system  is  to  pbtain  engine 
steering  and  shutoff  information  by  operating  on  measured  data.  These 
systems  are  basically  of  two  types,  radio  and  inertial.  A  summary  descrip 
tion  of  each  is  given  below. 

Radio  guidance  makes  position  and/or  velocity  measurements  from 
ground-based  radar.  These  measurements  are  input  to  a  ground-based 
computer  where  guidance  commands  are  computed.  These  commands  are 
sent  to  the  missile  control  system  via  a  radio  link.  A  radio  guidance  sys¬ 
tem  requires  missile  visibility  from  the  tracking  radar  site.  In  some 
cases  a  downrange  station  must  be  provided  to  insure  this  visibility.  A 
constraint  on  the  missile  attitude  may  also  be  imposed  because  the  beacon 
antenna  is  attached  to  the  missile.  The  accuracy  of  the  all -electronic 
radio  guidance  systems,  especially  those  with  long  base  lines,  is  generally 
better  than  may  be  obtained  with  inertial  guidance  systems  (IGS's). 

An  all-inertial  guidance  system  for  ballistic  missiles  makes  accelera- 
tion  measurements  in  a  gyro -stabilized  platform  coordinate  system. 

These  measurements  are  input  to  an  airborne  computer  where  guidance 
commands  are  computed  and  relayed  to  the  control  system.  The  control 
system  used  with  an  IGS  may  be  simplified  in  that  the  conventional  gyro 
attitude  reference  normally  used  with  radio  guidance  may  be  replaced  by 
the  inertial  platform.  The  electromechanical  sensors  used  in  inertial 
guidance  systems  are  subject  to  error.  The  precision  gyros,  used  in  the 
platform,  drift  over  long  periods  of  time  resulting  in  significant  coordinate 
system  errors.  The  bias  and  scale  factor  errors  in  the  accelerometer 
result  in  position  and  velocity  errors  in  the  navigation  loop.  The  above 
IGS  errors  are  offset  by  the  flexibility  offered  by  a  self-contained  non¬ 
radiating  system.  Present  day  IGS  accuracy  is  satisfactory  for  most 


Guidance  principles  that  are  directly  related  to  strapped-down  IG  systems 
will  not  be  discussed  in  this  report. 
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space  applications  in  the  earth's  gravitational  field;  this  system  may  also 
be  used  for  ballistic  missile  launch  guidance  to  lunar  impact.  For  more 
complicated  missions,  lasting  many  days  or  months,  an  IGS  may  be  aided 
with  a  star  and/or  planet  tracker.  Radio  guidance  may  also  be  used  as  an 
inertial  aid,  especially  when  target  acquisition  greatly  simplifies  the  guid¬ 
ance  task. 

The  principal  differences  between  radio  and  all-inertial  systems  then 
involve  the  form  of  the  original  measurements.  With  either  system  suf¬ 
ficient  data  must  be  measured  to  allow  the  calculation  of  position  and 
velocity  vectors,  R  and  R*  of  the  missile  in  some  suitable  coordinate 
system.  A  velocity  deviation  measured  from  the  desired  missile  velocity 
may  be  obtained.  The  thrust  acceleration,  may  then  be  steered  in 

direction  so  as  to  drive  the  yaw  and  pitch  components  of  the  velocity 
deviation  to  zero.  When  the  axial  (or  roll)  error  velocity  component  is 
reduced  to  zero,  engine  thrust  is  terminated.  The  missile  will  then 
coast  to  its  target  under  the  influence  of  mass  attraction  forces  exerted 
by  the  earth,  sun,  other  planets,  etc. 

Note  that  it  is  assumed  here  that  the  thrust  acceleration,  a.^,,  can 
only  be  steered  in  direction  and  not  controlled  in  magnitude.  This  lack  of 
control  does  not  impose  a  serious  restraint  because  the  magnitude  of  'a.p 
is  designed  to  be  reproducible  within  tolerances  of  a  few  percent. 

This  memo  will  be  concerned  with  all-inertial  guidance  systems  used 
for  ballistic  missile  guidance.  The  duration  of  powered  flight  in  a  ballistic 
trajectory  is  usually  several  minutes,  during  which  accelerations  may 
reach  15  g^'s.  All  guidance  equipment  on  board  the  missile  must  operate 
with  a  high  degree  of  reliability  in  this  environment. 

The  following  sections  of  this  niemo  will  discuss  inertial  coordinate 
systems,  measurement  errors,  the  navigation  loop,  guidance  computations, 
and  ballistic  missile  steering. 
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II.  INERTIAL  COORDINATE  SYSTEMS 

A  coordinate  system  (CS)  will  be  specified  in  this  section  in  terms  of 

s  \ 

three  mutually  orthogonal  inertial  directions.  These  directions  are  indi¬ 
cated  by  three  unit  vectors,  which  are  termed  basis  vectors.  Each 
coordinate  system  is  then  defined  by  specifying  the  directions  of  the  three 
basis  vectors.  Four  coordinate  systems  will  be  discussed.  The  orienta¬ 
tion  of  the  three  basis  vectors  will  in  each  case  be  apparent.  In  no  case 
will  the  transformation  from  one  CS  to  another  be  written  out  as  a  set  of 
equations.  The  orientation  of  the  basis  vectors  is  always  assumed  known, 
and  these  transformations  can  generally  be  written  by  inspection  in  terms 
of  the  geometry  of  the  system. 

The  inertial  measurement  unit  consists  of  three  accelerometers 
mounted  on  a  gyro -stabilized  platform.  In  this  memo,  it  will  be  assumed 
that  the  reference  platform  will  maintain  a  fixed  orientation  with  respect 
to  inertial  space.  Platform  hardware  such  as  gimbals,  the  porro  prism, 
etc.  ,  must  be  used  to  specify  three  orthogonal  platform  directions.  These 
directions  define  the  platform  coordinate  system.  Initial  platform  CS 
alignment  is  accomplished  with  respect  to  an  external  CS.  The  three  in¬ 
put  axes  of  the  platform  stabilizing  gyros  define  a  second  coordinate  sys¬ 
tem.  This  CS  must  be  defined  in  order  to  perform  platform  drift  error 
analyses.  It  will  be  shown  in  the  next  section  that  the  orientation  of  each 
gyro  must  be  specified.  The  three  accelerometers  are  fixed  rigidly  to 
the  inertial  platform.  The  directions  of  the  three  accelerometers  define 
the  accelerometer  coordinate  system.  This  is  an  important  CS  because 
it  is  the  one  in  which  the  measurements  of  sensed  acceleration,  'a.j,,  are 
made. 

A  fourth  coordinate  system  is  of  interest.  This  will  be  referred  to 
as  the  computing  coordinate  system,  and  is  the  one  used  for  guidance  cal¬ 
culations.  This  CS  may  be  an  earth -centered  inertial  CS  in  which  the 
expression  for  gravity  may  be  easily  calculated.  Other  coordinate  systems 
have  other  desirable  properties.  The  four  coordinate  systems. above  may  all 
be  equivalent  or  they  may  all  be  different  depending  on  actual  hardware  and 
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guidance  mechanizations.  In  any  case,  the  transformations  from  any  one 
CS  to  any  other  CS  will  be  known.  Hence,  for  purposes  of  this  memo,  it 
may  be  assumed  that  the  accelerometer  coordinate  system  and  the  com¬ 
puting  coordinate  system  are  identical. 

The  computing  CS  will  be  specified  in  terms  of  the  three  mutually 

■A  ^ 

orthogonal  unit  vectors  i  ,  j  ,  and  k  .  The  missile  position  and  velocity 

may  be  written  in  terms  of  this  set  of  inertial  basis  vectors  as 

^  -1 
R  =  X  i  +Y  j  +  Z  k 

— k  _1  -X 

R  =  Xi+Yj+Zk 


where  X  ,  Y,  Z  and  X,  Y,  Z  are  the  components  of  missile  position  and 
velocity  respectively.  The  computing  coordinate  system  inertial  basis 

^  -k 

vectors  i,  j,  k  may  be  used  to  evaluate  the  components  of  any  vector 
used  in  the  following  section.  For  example, 
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Ill,  INERTIAL  GUIDANCE  MEASUREMENT  ERRORS 


1 


Inertial  guidance  measurement  errors  affect  both  the  magnitude  and 
direction  of  the  sensed  acceleration.  These  error  sources  are  conven¬ 
iently  classed  as  follows: 

a)  Accelerometer  bias,  scale  factor,  and  second-degree  errors; 

b)  Accelerometer  and  platform  misalignments; 

c)  Fixed  and  acceleration-dependent  gyro  drift  rates. 

Errors  in  Class  (a)  affect  the  magnitude  of  the  sensed  acceleration, 
while  Class  (b)  and  (c)  errors  affect  the  direction. 

Detailed  discussions  of  inertial  guidai.-'e  error  analysis  techniques 
are  given  in  References  (1)  and  (2).  A  summary  description  is  given 
below. 

Sensed  Acceleration  Magnitude  Errors: 

All  accelerometers  contain  some  moving  parts.  When  acceleration 
is  sensed,  forces  within  the  instrument  cause  an  internal  displacement 
that  is  proportional  to  the  acceleration.  When  the  device  is  used  during 
ballistic  coast  periods,  no  such  forces  act  on  the  internal  parts  of  the 
instrument.  There  are  no  internal  displacements  in  a  zero-g^  field  when 
changes  in  "g  within  the  instrument  are  neglected.  Hence,  the  reading  of 
an  ideal  accelerometer  is 


! 


where 

lA 

R 


_\ 

g 


acceleration  of  the  missile  with  respect  to  inertial  space 
acceleration  of  gravity 

the  sum  of  all  external  accelerations  acting  on  the  missile 
except  g 
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Let  and  denote  the  ideal  and  actual  readings  of  the  i 

accelerometer.  The  accelerometer  error  Aa_,  is  then 

Tx 


Aa 


Tx 


^Tx  ■  ^Tx 


Accelerometer  measurement  errors  result  from  bias,  cross  acceleration, 
and  nonlinearity  errors.  A  general  expression  for  accelerometer  error 
is  given  as 


^^Tx  ■  ^o  ^1  ^Tx  ^2  ^Ty  ^3  ^Tz  ^1 1  ^Tx^  ^12  ^Tx^’Ty 

(3.2) 

where 


*"1 
k^j  k^ 

^11 
k etc. 


^Tx’^Ty’^Tz 


bias  coefficient 

scale  factor  coefficient 

cross  axis  coefficients 

second -degree  error  coefficient 

additional  nonlinearity  error  coefficients 

components  of  sensed  acceleration 


Expressions  similar  to  equation  (3.  2)  are  used  for  the  j  and  k  accelero¬ 
meters.  The  relative  magnitude  of  the  error  coefficients  in  equation 
(3.  2)  depends  on  the  particular  accelerometer  design.  The  actual  values 
are  usually  obtained  from  the  manufacturer  or  impartial  testing  agencies. 


Accelerometer  and  Platform  Misalignments: 

The  accelerometers  cannot  be  mounted  on  the  platform  without  intro¬ 
ducing  some  error.  Similarly  the  platform  cannot  be  leveled  and  aligned 
in  azimuth  without  error.  As  a  result,  six  misalignment  angles,  two  for 
each  accelerometer,  are  required  to  express  the  misalignment  of  the 
accelerometer  package  with  respect  to  inertial  space. 
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Let  i,  j,  k  denote  the  desired  orientation  of  the  three  accelerometers. 
The  actual  orientations  i',  j',  and  k'  are  then 

^  =  "i  +  <t>i^  +  013^ 

?  =  T  +  ^23^ 

k'  =  k  +  +  032T 


where  the  0^2’  ‘  ’  ^32  small  misalignment  angles.  For 

example,  0^3  misalignment  of  the  accelerometer  in  the  k 

direction.  The  misalignment  angles  may  be  specified  in  several  ways. 

For  example,  three  angles  may  be  used  to  specify  the  platform  level  and 
azimuth  alignment.  These  angles  are  equivalent  to  ^^31  ^33*  ^12’ 

The  i  accelerometer  may  then  be  mounted  on  the  platform  with  zero 
error  assumed.  The  j  accelerometer  is  then  mounted  on  the  platform 
with  respect  to  i;  the  error  is  02i‘  ^  accelerometer  is  then 

mounted  with  errors  ^3^  and  ^32*  Again,  let  a,j,y, ,  and  a,p^ 

denote  the  ideal  accelerometer  measurements.  The  actual  measurements 
^Tx’  ^Ty’  ^Tz’  obtained  from  equation  (3.3)  as 


'Tx 


^Tx  ^12  ^Ty  ^ 


13  ^Tz 


^Ty  ■  ^Ty  ^ 


23®^Tz 


+  0 


21  ^Tx 


(3.4) 


^Tz  ^Tz  "^31  ^Tx  ^32  ^Ty 

The  last  two  terms  in  each  of  equations  (3.4)  are  error  terms.  These 
are  easily  calculated  using  standard  missile  simulation  programs. 


Platform  Drift  Error: 

Coordinate  system  errors  result  from  drifts  in  the  platform  stabi¬ 
lizing  gyros.  In  addition  to  the  steady  drift  rate  of  the  gyro,  acceleration- 
dependent  drift  rates  may  exist.  These  latter  drift  rates  are  referred  to 
as  mass  unbalance  and  anisoelastic  drift  rates.  The  following  discussion 
is  included  to  aid  in  understanding  the  nature  of  IGS  platform  drift. 
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A  platform  uses  three  single-degree-of-freedom  gyros  or  two  two- 
degree -of-freedom  gyros.  The  three  orthogonal  axes  of  a  single-degree- 
of-freedom  gyro  are  shown  in  Figure  3,  1, 


©^(output  axis) 


Figure  3.  1  -  Axes  of  a  SDF  Gyro 

In  Figure  3.  1,  the  subscript  n  =  1,  2, or  3  is  added  to  S,  I,  and  O 

th 

to  indicate  the  spin,  input,  or  output  axis  of  the  n  platform  gyro,  A 
platform  SDF  rate  integrating  gyro  is  designed  to  measure  only  angular 
rates  about  the  input  axis.  An  angle  pickoff  on  the  output  axis  measures 
the  input  angular  rate.  Angular  rates  about  the  output  axis  do  not  cause 
a  pickoff  indication  because  the  viscous  fluid  damping  constrains  the  case 
motion  to  follow  the  inner  gimbal  motion. 

A  conventional  two-degree-of-freedom  gyro  differs  from  a  SDF  gyro 
in  that  an  additional  gimbal  and  pickoff  are  added.  As  such,  the  2-DF  gyro 
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may  be  considered  equivalent  in  performance  to  two  SDF  gyros.  When  two 
2-DF  gyros  are  used  to  stabilize  the  platform  a  redundant  control  direction 
is  provided.  This  redundant  direction  will  be  ignored  in  this  analysis.  The 
subscript  n  will  be  used  to  denote  any  one  of  three  orthogonal  input  axes 
that  are  used. 

A  2-DF  gyro  does  not  require  the  viscous  damping  restraint.  In  fact,  the 
flotation  fluid  has  negligible  damping.  The  additional  gimbal  and'pickoff 
is  provided  so  that  angular  displacements  about  both  input  axes  may  be  meas¬ 
ured.  The  directions  of  these  axes  are  shown  in  Figure  3.2.  Note  that 
each  input  axis  is  also  an  output  axis. 


I 


2 


Figure  3.2  -  Axes  for  a  2-DF  Gyro 


The  angular  momentum  of  the  gyro  spin  rotor  is  large  compared  to 
the  gimbal  momenta  and  certain  other  dynamical  terms  discussed  in 


/ 


References  (3)  and  (4).  Hence,  the  law  of  the  gyro  is  approximated  as 


T  =  w  A  H 
n  n  n 


(3.5) 


where 


=  vector  angular  momentum  of  the  gyro  wheel 
=  residual  unbalance  torques  within  the  gyro 

The  symbol  A  is  used  to  denote  the  vector  cross  product. 

Both  SDF  and  2-DF  gyros  drift  as  a  result  of  unbalance  torques  acting 
about  the  output  axis.  When  both  sides  of  equation  (3.  5)  are  multiplied  by 
(scalar  product)  the  result  is 


T  •  O 
n  n 


(co  aH  )  •  O 
n  n  n 


H  u  •  (S  AO  ) 
n  n  n  n 


=  -H  w  .1 
n  n  _i 


Hence,  the  platform  drift  rates  about  the  three  input  axes,  u  ,  may  be 
sl<  ^ 

written  as 


^  T  .  O 
H  n  n 
n 


(3.6) 


Fixed,  i.  e.  ,  constant  or  steady,  gyro  drift  rates  are  usually  denoted 

by  E..  These  drift  rates  result  from  fixed  residual  torques,  T^,  within 

the  gyro  such  as  bearing  frictions  and  pickoff  reactions.  Then  the  fixed 
tVi 

drift  rate  of  the  n  platform  gyro  is 


R  =  -  7^  T  .  O 
n  H  on  n 

n 


(3.7) 


The  presentation  in  this  section  is  based  on  original  work  by  T.  W.  Layton 
and  H.  Cohen. 
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i‘Vk 

A  mass  unbalance  within  the  n  gyro  is  approximately  equivalent  to  the 
torque 


T  =  m  6  A  arri 
n  n  n  T 


(3.8) 


whe  r  e 

m^  =  internal  gyro  mass 

6^  =  vector  displacement  of  the  center  of  mass  of  the  rotor  gimbal 
system 

a,p  =  sensed  acceleration  vector 

The  quantities  6^  and  a^  may  be  written  as 

"t  =  °n  +  4  •  y  X.  "  |"t  •  ®n 

^  ^  -X 

6  =  (6  .0)0  +  (6  •  I  )  I  +  (6  .  S' )  S 

n  '  n  n'  n  '  n  n  n  '  n  n'  n 

Hence,  the  n^^^  gyro  drift  rate  resulting  from  mass  unbalance  is 


=  U  (a^  .  SJ  -  Ug  (a^  • 
n  ^  n 


(3.9) 


where 


m 


U, 


n 


n 


n 

H 

n 

(5 

n 

m 

n 

K 

•t> 

H 

n 

(3.9A) 

(3.9B) 


Anisoelastic  gyro  drift  rate  results  when  the  gyro  is  placed  in  an 

acceleration  field.  In  this  environment  the  gyro  will  deform.  The  result- 

ing  mass  shift  is  obtained  as  m  €  ,  where 
®  n  n 

-  Ko„<"t  •  °n  "  •  y  t  +  ■'sn  <"t  '  =«>  10) 
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th 

The  constants  ^In'  coefficients  of  the  n  plat¬ 

form  gyro.  The  resulting  residual  torque  is 


n  n  n  T 


(3.11) 


Substituting  equations  (3.  11)  and  (3.  10)  into  (3.6)  gives 


where 


“A  = 

n_ 


2K 

n 


n  n 


H 

n 


(3.12) 


(3.  12A) 


K 

n 


anisoelastic  drift  rate  coefficient 
for  the  n^^  platform  gyro 


The  gyro  drift  rate,  “An’  referred  to  as  the  anisoelastic  drift  rate  of 
the  n^^  platform  stabilization  gyro.  The  gyro  is  said  to  be  isoelastic  if 
This  gyro  design  goal  has  not  yet  been  achieved. 

The  resultant  gyro  drift  rate  is  denoted  by  iTj^.  An  expression  for 

(r_  is  obtained  by  adding  the  fixed,  mass  unbalance  and  anisoelastic  drift 

til 

rates.  The  result,  for  the  n  platform  gyro,  is 


o-_  =  R  +  U,  a^  -  U,  a,  +  2K  a^  ar.  (3.13) 

D  nIS  II  nIS 

n  n  n  n  n  n 

where 


ag  a.j, 

.  “s 

n 

(3.13A) 

n 

(3.  13B) 

af  =  a-rj. 

n 

•  I 

n 

The  above  discussion  applies  to  conventional  ball-bearing  floated  plat¬ 
form  gyros.  Gyros  utilizing  other  designs  may  have  internal  acceleration- 
dependent  torques  that  are  not  calculated  above. 
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There  are  many  ways  that  three  SDF  gyros  may  be  oriented  on  the 
platform.  One  such  orientation  is  given  in  Table  3.1,  where  it  is  assumed 
that  the  thrust  acceleration  is  always  contained  in  the  X-Z  plane. 


Gyro  1 

Gyro  2 

Gyro  3 

X  Direction 

1^2  cos  a  +  XS^  sin  a 

°3 

Y  Direction 

s, 

I, 

s. 

1 

2 

3 

Z  Direction 

-O, 

-St  sin  a  +  0_  cos  a 

I- 

1 

2  2 

3 

Table  3.  1  -  An  Optimum  Gyro  Orientation 

The  following  features  of  Table  3.  1  should  be  noted: 

a)  The  input  axes  are  orthogonal  and  define  the  gyro  coordinate 
system. 

b)  The  SDF  gyros  1  and  3  could  be  replaced  by  a  single  2-DF  gyro. 

c)  Since  a^  is  in  the  X-Z  plane,  there  is  no  anisoelastic  drift  rate 
regardless  of  the  value  of  the  angle  a. 

d)  The  platform  could  be  rotated  about  the  j  direction  through  an 
angle  p  without  introducing  anisoelastic  drift  rates. 

The  gyro  orientation  given  in  Table  3.  1  successfully  eliminates 
anisoelastic  drift  rates,  but  not  mass  unbalance  drift  rates.  It  is  possible 
to  choose  the  ngles  a  and  p  to  minimize  these  latter  drift  rates.  When 
this  has  been  done  the  gyro  orientation  is  said  to  be  optimized. 
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IV.  POSITION  AND  VELOCITY  DETERMINATION 


A  ballistic  missile  guidance  system  does  not  always  require  explicit 
information  as  to  position  and  velocity,  as  in  the  case  of  the  Q-guidance 
method  proposed  by  MIT.  The  guidance  methods  to  be  discussed  in  this 
memo,  however,  require  that  the  inertial  position  vector,  R  ,  and  the 
inertial  velocity  vector,  R  ,  be  determined.  The  equation  to  be  solved 
is 

R  (t)  =  g  (R)  +  Sirj,  (t)  (4.  1) 

where 


g  =  gravitational  acceleration 

a.rj,  =  sensed  acceleration  (measured  by  the 
accelerometers) 

R  =  inertial  acceleration 


Equation  (4. 1)  is  integrated  in  the  IGS  computer.  The  block  diagram 
shown  in  Figure  4.  1  indicates  the  steps  in  the  solution.  The  following 
notation  is 


Figure  4,1-  Position  and  Velocity  Determination 
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used  in  Figure  4.  1 

=  initial  position  vector  of  missile  ^ 

=  initial  velocity  vector  of  missile 

Equation  (4.  1)  is  solved  in  the  computing  coordinate  system.  If  the 
accelerometer  coordinate  system  is  not  the  same  CS,  a  transformation  of 
aj  to  the  computing  CS  must  be  included. 

Division  and  square  foot  operations  are  slow  processes  for  many 
IGS  computers.  Hence,  the  direct  computation  of  g  is  a  slow  process 
when  "g  has  the  form 

t  =  -  ™  «  (4.2) 

r 

where 

R  =  missile  position  vector  measured  from  the 
center  of  the  earth 


In  equation  (4.2)  the  position  vector  of  the  missile,  R,  is  expressed 
in  a  coordinate  system  with  origin  at  the  center  of  the  earth.  At  worst 
this  requires  a  simple  transformation  to  the  computing  CS.  The  gravity 
function  may  be  expanded  into  a  Taylor's  series  of  the  form 

g  =  -  (R  +  RJ  F  (X,  Y,  Z)  (4.3) 

where 

F  =  C  V  C,  X  +  C,  Z  +  C.X^  +  C.Z^  +  C.XZ 
o  1  2  4  5  o 

A  complete  description  of  this  expansion  is  given  in  Reference  (5). 
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In  equation  (4.  3)  it  is  assumed  that  the  nominal  thrust  acceleration 
vector  is  always  contained  in  the  X  -  Z  plane.  The  above  computational 
form  for  g  includes  oblateness  terms  while  the  form  given  in  equation 
(4.2)  does  not.  The  coefficients,  C^,  required  by  equation  (4.3)  are 
obtained  by  Taylor's  series  expansion  about  a  downrange  expansion  point. 
Other  methods  for  mechanizing  the  gravity  computation  exist.  For  example, 
equation  (4.2)  may  be  calculated  using  difference  equations.  In  some  cases, 
for  a  short  high-altitude  burn  it  will  suffice  to  use  a  constant  vector  for 
gravitational  acceleration. 

Most  of  the  accelerometers  currently  used  in  IG  systems  are  integrat¬ 
ing  accelerometers.  That  is,  the  output  of  the  accelerometer  is  thrust 
velocity,  V.^,,  where 


When  this  is  the  case  the  navigation  loop  shown  in  Figure  4.  2  may  be  used. 


Figure  4.2  -  A  Possible  Navigation  Loop  When  an  Integrating 
Accelerometer  Is  Used. 
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The  initial  position  and  velocity  and  required  in  Figures  4.  1 

and  4.2  are  readily  determined  when  the  missile  is  launched  from  the 
earth.  In  some  space  missions  a  second  burning  peripd  may  be  required 
after  the  coast  (ballistic  free  flight)  phase.  The  determination  of  initial 
conditions  for  this  second  burn  may  be  made  in  two  ways.  Since  ‘g  is  the 
only  external  acceleration  acting  on  the  missile  during  coast,  it  is  only 
necessary  to  integrate  "g  (twice)  during  the  coast  period.  When  this 
method  is  used  the  accelerometer  inputs  in  Figures  4.  1  and  4.2  are  opened 
so  that  accelerometer  bias  errors  are  not  introduced. 


The  second  method  requires  that  the  six  initial  conditions  at  the  start 

of  the  second  burn  be  predicted  using  first-burnout  position  and  velocity 

information.  It  is  also  necessary  to  predict  the  time  of  free  flight,  t^^, 

in  order  to  provide  an  engine  start  command  for  the  second  burn.  It  is 

shown  in  Section  V  that  at  first  burnout,  the  missile  will  satisfy  four 

guidance  restraints  at  the  start-point  of  the  second  burn.  This  start- 

point  is  the  "target"  for  the  first  burn.  When  the  four  guidance  restraints 

•  « 

are  X  =  X.p,  Y  =  Z  =  Zj,  Z  =  Zj,  then  three  quantities  remain 

to  be  determined.  These  may  be  predicted  using  expressions  of  the  follow¬ 
ing  form: 


where 


^T 

Ci^AY 

+ 

CjjAZ 

+  Cj^At 

""  ^TN 

C2,AX  + 

C^^AY 

+ 

+  C24'“ 

•if  =  ‘ffn  *  •=41'^^  +  <^42^^  *  ^43^^  *  ‘=44'“ 

fix  =  X  -  Xj, 

=  Y  - 

AZ  =  Z  -  z„ 

N 

At  =  t  .  tj, 
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^N’  '^N*  ^N’  ^ 


nominal  position  and  time  at  first  burnout 


X.  Y.  Z.  t 

X  Y 
TN’  TN 


actual  position  and  time  at  first  burnout 
nominal  velocity  components  at  the  target 
nominal  time  of  free  flight 


When  the  above  method  is  used  the  guidance  computer  (except  the 
clock)  as  well  as  the  accelerometer  and  accelerometer  electronics  may  be 
shut  off.  This  results  in  a  saving  of  power,  and  therefore  weight. 
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V.  THE  REQUIRED  VELOCITY  VECTOR 

A  ballistic  missile  engine  whose  thrust  may  be  controlled  in  magnitude 
as  well  as  direction  is  said  to  have  thrust  magnitude  control.  Ballistic 
missiles  with  this  feature  can  be  made  to  fly  the  nominal  desired  time- 
specified  trajectory.  That  is,  such  missiles  could  be  controlled  to  burn¬ 
out  at  a  prespecified  position,  velocity,  and  time. 

It  will  be  assumed  in  this  memo  that  the  ballistic  missile  engines  used 
do  not  have  thrust  magnitude  control.  However,  it  will  be  assumed  that 
the  thrust  and  mass  magnitudes  are  reproducible  within  a  few  percent. 
These  ballistic  missiles  can  then  be  steered  only  by  pitching  and  yawing 
the  thrust  vector.  With  this  type  of  control  it  is  possible  to  steer  the 
missile  velocity  vector  to  the  value  specified  by  some  function  of  missile 
position,  R,  and  time,  t.  This  desired  velocity  function  will  be  termed 
the  required  velocity  vector,  Vj^  =  Vp^  (R,  t).  This  vector  function  is 
obtained  on  the  limiting  assumption  that  an  impulsive  thrust  engine  is 
used,  That  is,  at  any  time  during  the  powered  flight,  it  is  assumed  that 
a  velocity  increment,  AV,  may  be  added  such  that 


k  -A 

+  AV  = 


This  is  referred  to  as  a  limiting  assumption  because  additional  guidance 
constraints  may  be  imposed  when  the  actual  finite  thrusting  capabilities 
of  the  missile  are  used  in  the  analysis. 

The  limiting  assumption  will  be  made  in  the  analysis  that  follows. 
Then  the  specification  of  four  general  guidance  constraints  is  necessary 
in  order  to  uniquely  determine  Vp  (R,  t).  A  consideration  of  Figure  5.  1 
will  help  to  clarify  this  matter. 
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Figure  5. 1  -  Three  Coast  Trajectories,  All  Hitting 
the  Desired  Target 


The  three  free -flight  trajectories  shown  in  Figure  5,  1  all  "hit”  the  de¬ 
sired  target.  That  is,  after  a  period  of  time,  t^^,  each  of  the  three 
trajectories  satisfies  the  three  guidance  constraints  that 


Consider  the  two  trajectories  originating  at  point  A.  In  order  to  specify 
the  required  velocity  vector  at  point  A,  an  additional  guidance  constraint 
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must  be  imposed.  For  example,  it  may  be  required  to  hold  the  time  of 
free  flight  to  a  constant  value,  as  is  done  in  trajectories  A  and  B. 
Equally  well,  any  one  of  the  following  restraints  could  be  prespecified: 

a)  Burnout  velocity  magnitude 

b)  Burnout  vertical  velocity 

c)  Burnout  velocity  vector  elevation  angle 

d)  Burnout  energy 

e)  Burnout  angular  momentum  magnitude 

f)  Velocity  magnitude  at  the  target 

g)  Vertical  velocity  at  the  target 

h)  Total  time  of  flight  from  liftoff  to  the  target 

In  addition  to  the  above,  many  more  restraints  which  are  desirable 
could  be  listed. 

It  is  not  necessary  to  include  hitting  the  target  as  three  restraints. 

a 

A  required  velocity  vector  field,  Vj^  (R.*t),  can  be  found  that  will  cause 
the  free-flight  velocity  vector  to  satisfy  the  following  four  conditions; 

1)  X  =  X,j, 

Z)  Z  = 

3)  Y  = 

4)  Constant  total  time  of  flight  from  liftoff 

•  • 

The  t^j,  X,j,,  Zrj,,  and  Y,p  will  all  vary  in  the  preceding  example 

depending  on  the  burnout  time  and  position  vector.  Another  interesting 
example  will  be  included.  It  is  possible  to  steer  the  missile  so  that  the 
velocity  vector  at  burnout  has  some  desired  value.  In  this  case  the  four 
guidance  restraints  are  the  following: 
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4)  t. 


In  this  example  the  burnout  position  and  time  will  vary  depending  on  the 
off  nominal  performance  of  the  propulsion  system,  etc.  The  above  ex¬ 
amples  serve  to  illustrate  that  all  essential  guidance  information  is 
included  in  the  specification  of  the  required  velocity  vector.  .  This 

follows  from  the  fact  that  all  missiles  can  be  controlled  (theoretically) 
such  that  they  will  have  exactly  the  required  velocity  vector  at  burnout. 
They  will  therefore  satisfy  the  prescribed  guidance  restraints.  Any  two 
missiles  that  burn  out  at  the  same  vector  position  and  time  will  have 
identical  free -flight  trajectories.  Note  that  the  missile  must  be  steered 
in  pitch  in  order  to  achieve  (R.t), 

The  required  velocity  vector  is  sometimes  discussed  with  a  meaning 
different  from  that  described  above.  For  example,  the  four  guidance 
restraints  may  be  the  following: 

1,  2.  3)  R  = 


4)  Arbitrary  vertical  velocity,  Z,  at  burnout. 


Note  that  the  fourth  restraint,  above,  is  quite  different  from  condition 
b)  previously  discussed,  where  it  was  required  that  the  burnout  vertical 
velocity  be  steered  to  a  prespecified  constant.  Condition  4)  is  used 
whenever  it  is  required  to  pitch  the  missile  according  to  some  preassigned 
pitch  program.  In  this  case  the  missile  vertical  velocity,  Z,  will  vary 
according  to  missile  performance.  The  variation  in  Z  maybe  indicated 
by  writing  the  guidance  control  function  as: 
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Vr  =  Vr  (R.  Z,  t) 


where  Z  is  arbitrary. 


Using  Vr  in  this  form  leads  to  satisfactory  performance  for  missiles 
that  are  designed  to  impact  on  the  earth.  Note  that  only  three  guidance 

-i 

restraints  may  now  be  imposed.  That  is,  at  R  =  R,p  along  the  coast 
trajectory,  the  coast  time  of  free  flight  will  necessarily  vary  with 
missile  performance.  This  memo  will  not  further  discuss  guidance 
systems  using  Vr  =  Vr  (R,  Z,  t).  Only  guidance  systems  using  Vr 
=  Vr  (R,  t),  which  require  pitch  steering,  will  be  discussed. 

Some  simple  examples  dealing  with  required  velocity  vector  calcula¬ 
tion  will  now  be  given.  The  calculation  of  V-  involves  only  the  free- 
flight  equations  of  motion.  On  a  flat  earth  these  may  be  written  as 


R  =  g 

where  g  is  a  constant.  The  solution  is 


(5.  1) 


R  =  R  +  R  t  +  4-  gt^ 
o  o  2  ^ 


(5.2) 


In  equation  (5,  2),  R  and  R  denote  the  burnout  conditions.  It  is  now 

o  o 

required  to  find  the  burnout  velocity  that  will  cause  the  missile  to  coast 
to  the  target,  R,^,  (t),  after  a  time  of  free  flight,  t^^.  Since  it  is  assumed 
that  the  missile  can  burn  out  at  any  position,  R,  equation  (5,2)  can  be 
written  as 


R,p(t)  =  ^  ^  J  S  ^ 


if 


(5.2A) 


The  target  position  is  written  as  R,^  (t)  to  indicate  that  the  position  of 
the  desired  target  may  vary  with  time,  The  solution  for  V_  on  the  flat 
earth  is  then 


VR(R,t)  = 


RT(t)  -  R 


-Tgt 


ff 


(5,  2B) 
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If  the  t££  is  now  chosen  as  the  fourth  guidance  restraint,  then  this  pre¬ 
specified  value  (a  constant)  may  be  substituted  in  equation  (5.2B)  resulting 
in  the  desired  expression  for  Vj^.  If  the  fourth  guidance  restraint  is  that 
the  total  time  of  flight  should  be  a  prespecified  constant,  T,  then 
tff  =  T  -  t  should  be  substituted  into  equation  (5. ZB), 

It  may  be  required  that  the  vertical  velocity  of  the  missile  at  the 

target  be  some  pre specified  value,  Z_.  Then,  when  g  =  g  c  0, 

1  X  y 


,:g  -  32 


(5.3) 


When  is  eliminated  from  equation  (5.  3)  and  the  Z-component  of 

equation  (5.  2B),  an  expression  for  time  of  free  flight  is  obtained  as 


Z^^2g{Z,j.  -  Z) 


(5.4) 


The  desired  expression  for  required  velocity  is  now  obtained  by  substituting 
t££  from  equation  (5. 4)  into  (5.  2B), 

In  realistic  cases  must  be  obtained  for  oblate  spheroids.  This 

solution  can  be  obtained  with  sufficient  accuracy  by  analytical  methods. 
However,  these  solutions  will  involve  square  roots  and  divisions  which 
are  difficult  operations  for  many  IGS  computers.  Experience  indicates 
that  a  loss  in  accuracy  results  when  these  analytical  solutions  are  expanded 
in  power  series  expansions.  Hence,  curve  fitting  techniques  (least  squares) 
are  employed  using  high-speed  digital  computers.  It  is  desired  to  find 

-Jk  -X 

Vj^(R,t)  in  the  following  form: 

=  Sqq  +  SjqAX  +  s^qAY  +  a^gAZ 

+  a^Q  At  +  ajjAX^  +  AX  AY 


+  . . .  +  a^^  At 


(5.5) 


-24- 


There  are  15  terms 


Similar  expansions  are  obtained  for  and 

in  each  expansion.  In  equation  (5.5), 


AX  =  X  -  X 


N 


AY  =  Y  -  Y 


N 


AZ  =  Z - Z 


N 


At  =  t  -  t 


N 


Xj^,  Yj^,  Zj^,  tj^  =  nominal  (standard)  burnout  position  and  time 
a^j  =  guidance  constants,  to  be  determined 
The  a^j  in  equation  (5.5)  are  analogous  to  partial  derivatives,  as 


8X 


^10  -  W 


Burnout 


= 

12  9X  9Y  /  Burnout 


except  that  these  "partial  derivatives"  are  obtained  by  least  squai’es 
techniques.  Note  that  apQ  in  equation  (5.5)  is  the  nominal 
value  of  the  X- component  of  the  burnout  velocity  vector  or 


^00  ^RxN 


No  third-order  terms  are  given  in  equation  (5.  5).  In  general,  the 
listed  terms  have  been  found  to  result  in  negligible  fitting  errors  over 
the  expected  burnout  dispersion  box. 

The  burnout  dispersion  box  is  defined  by  the  expected  variations  in 
AX,  AY,  AZ,  and  At  at  burnout.  These  variations  result  from  varia¬ 
tions  in  the  thrust,  specific  impulse,  initial  mass,  etc. 
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Figure  5.  2  -  Nominal  and  Nonstandard  Required 
Velocity  Vectors  at  Burnout. 


The  dispersion  box  is  shown  in  Figure  5,2  along  with  nominal  and 
nonstandard  required  velocity  vectors  at  burnout.  The  size  of  the  _ 
burnout  dispersion  box  depends  principally  on  the  propulsion  and  control 
system  tolerances.  The  shape  depends  on  the  guidance  restraints  chosen 
and  the  type  of  missile  steering  employed  by  the  IGS, 

In  order  to  determine  the  a^j  in  equation  (5,  5)  by  least  squares  pro¬ 
cedures,  it  is  necessary  to  have  input  data.  To  obtain  these  data,  it  is 
first  necessary  to  estimate  the  size  and  shape  of  the  dispersion  box. 

When  this  has  been  done  a  large  number  of  points  may  be  designated  with¬ 
in  the  dispersion  box  either  at  random  or  by  some  systematic  procedure. 
Any  one  of  many  iterative  techniques  may  now  be  used  to  determine  the 
required  velocity  vectors  at  the  designated  points.  With  some  techniques, 
as  many  as  15  machine  runs  may  be  required  to  determine  Vj^  at  one 


1 

! 
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point.  By  definition,  the  resulting  vector  field,  Vj^  (R,  t)  must  satisfy 
the  four  chosen  guidance  constraints.  This  vector  field  is  then  the  input 
data  to  the  least  squares  routine. 

A  second  method  for  obtaining  these  data  has  been  used  with  success 
at  STL.  This  method  has  the  advantage  that  iterative  techniques  are  not 
required.  Assume  that  the  four  guidance  restraints  are 

1.2.  3)  R  =  Rr^ 

4)X  =  Xrj, 

s  e 

Then  and  Z,p  (as  well  as  the  t,,)  will  vary.  The  free  flight  may 

•  • 

be  started  at  the  target  position  with  X  =  X„.  Runs  may  now  be  made, 
integrating  backwards.  These  runs  require  that  Y,j,  and  Z,p  be 
varied  in  amounts  depending  on  the  size  of  the  expected  dispersion  box. 
When, the  missile  has  coasted  (backwards)  for  the  nominal  time  of  free 
flight,  the  quantities  R,  R,  and  t££^  are  recorded.  These  data  are  also 
recorded,  along  the  same  free-flight  trajectory  at  the  times 

t  =  t££j^  ±  n  At  (n  =  1,2,  .  .  .  ,  N) 

where  At  is  also  dependent  on  the  size  of  the  dispersion  box.  If  N  =  3, 

then  seven  sets  of  data,  i.  e.  ,  seven  required  velocity  vectors,  are  obtained 

•  • 

for  each  run.  Then  from  the  nominal  and  two  variations  on  Y>j,  and  Zrj,, 

35  sets  of  data  will  be  obtained,  which  may  be  input  into  the  least  squares 
routine.  Note  that  it  is  desirable  to  define  all  guidance  restraints  at  the 
target,  when  using  this  method.  If  the  velocity  magnitude  at  the  target 
were  constrained,  then  X.^,  could  be  solved  for  and  the  above  discussion 

applies.  If  the  total  time  of  flight,  T,  were  constrained  (as  condition  4) 

*  •  • 

it  would  be  necessary  to  vary  AX.^,,  as  well  as  AY,j.  and  AZ.^..  The 
quantity  At^^,  appearing  in  the  expansion  would  then  be  eliminated  using 
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or 

-  At 

~  ^*ff  ■  ^ff  ■  *ffN 

where 

At 

= 

Methods  for  steering  the  missile  such  that  R  will  be  controlled  to 

be  equal  to  Vn  (R>  t)  will  be  presented  in  the  next  section.  Note  that 

such  control  is  equivalent  to  controlling  V  (R,  R,  t)  to  zero, 

S 

where 


V 

g 


(5.6) 


The  quantity  V  is  used  as  tKe  basic  guidance  control  function  in  the 
g 

MIT  guidance  system. 
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VI.  BALUSTIC  MISSILE  STEERING  METHODS 

A  ballistic  missile  is  steered  by  controlling  the  direction  of  the 
thrust  vector.  For  flight  through  the  atmosphere,  the  missile  must  be 
steered  so  that  a  gravity  turn  is  approximated.  A  gravity  turn  results 
when  the  missile  angle  of  attack,  and  hence  aerodynamic  forces  normal 
to  the  missile  axis,  vanish.  This  requirement  is  made  because  normal 
forces  of  this  nature  set  up  bending  moments  along  the  missile  axis  that 
could  cause  the  missile  to  break.  It  is  shown  in  Appendix  B  that  a 
gravity  turn  will  result  when  §  is  chosen  as 

5  R  -  CO  j,  A  R  +  W 

where  R  is  the  missile  inertial  velocity,  coj,  is  the  earth  angular 
velocity,  and  W  is  the  wind  velocity.  The  wind  velocity  is  usually 
neglected. 

It  will  be  assumed  in  this  section  that  the  commanded  attitude  turning 
rates  required  for  the  gravity  turn  are  obtained  from  a  simulation  incor¬ 
porating  the  above  equation.  The  required  turning  rates  may  then  be 
printed  out  and  stored  in  a  missile  programmer.  Alternate  methods  for 
accomplishing  the  gravity  turn  have  been  studied.  For  example,  the 
required  attitude  may  be  expressed  a§  an  empirical  function  of  actual 
missile  position  and  velocity.  This  procedure  minimizes  the  dispersions 
that  result  v/hen  closed-loop  guidance  steering  is  initiated. 

The  remaining  paragraphs  in  this  section  will  assume  that  the  gravity 
turn  has  been  completed  and  that  the  missile  is  out  of  the  atmosphere. 

Hence,  no  aerodynamic  forces  will  be  considered  to  act  on  the  missile. 

Ballistic  missiles  generally  have  long,  slender  configurations  with  the 
thrust  chamber  gimbal  located  many  feet  away  from  the  center  of  gravity 
(C.G.  ).  Since  the  consideration  of  missile  control  systems  would  tend  to 
mask  the  guidance  problem,  the  following  discussion  will  be  limited  to 
point  mass  missiles.  Such  missiles  are  approximated  by  cylindrical 

The  symbol  ~  denotes  parallelism. 
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shapes  with  the  thrust  chamber  gimbal  located  at  the  C.G,  The  direction 

of  the  thrust  attitude  vector  will  be  denoted  by  |  in  this  section.  The 

IGS  is  also  assumed  to  be  located  at  the  C.  G.  The  purpose  of  the 

IGS  is  to  determine  a  commanded  value  of  the  thrust  attitude  vector, 

-u  . 

such  that  the  missile  velocity  will  be  controlled  to  the  required 
velocity  vector. 

Steering  methods  will  be  developed  in  this  section  for  point  mass 
missiles.  These  methods  will  be  extended  in  Section  VII  to  actual  missiles, 
where  the  differential  equations  relating  ^  and  will  be  obtained. 

It  is  sufficient  for  the  present  purpose  to  note  that  ^  will  lag  as  a 
result  of  control  system  reaction  time,  engine  gimbal  Inertia,  and  any 
limiting  that  might  be  present.  It  is  assumed  here  that  the  control  sys¬ 
tem  response  is  rapid  and  that  these  lags  may  be  neglected.  However, 
thrust  misalignments  are  considered  to  exist,  and  integral  control  terms 

are  added  to  reduce  their  effect.  Hence,  ^  and  |  will  be  treated  as 

c 

identical  quantities  in  this  section. 

The  problem  of  steering  the  missile  may  be  approached  in  many 
ways.  The  solution  to  this  problem  is  not  unique.  The  methods  pre¬ 
sented  in  this  section  have  been  found  to  be  satisfactory,  but  no  claim  is 
made  that  these  methods  are  optimum.  A  short  discussion  of  the  nominal 
trajectory,  follows. 

Trajectories  for  nominal  missiles  are  usually  generated  by  choosing 
a  thrust  attitude,  ^j^(t),  such  that 

a)  R  (R,  t) 

b)  Minimum  propellant  is  wasted 

c)  Heating  constraints,  etc. ,  are  satisfied 

d)  Gravity  turn  is  flown  through  the  atmosphere 

It  may  require  many  runs  on  a  large  digital  computer  of  the  powered- 
flight  simulation  to  determine  subject  to  the  listed  restraints. 
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The  nominal  thrust  attitude  is  determined  from  liftoff  in  the  powered- 

flight  simulation.  Note  that  will  change  in  direction  except  for  verti- 

^  i 

cal  flight.  It  is  customary  to  define  a  constant  vector  Pj^  such  that 

The  plane  defined  by  is  said  to  be  the  pitch  plane. 

Since  guided  flight  in  a  vacuum  will  be  assumed,  only  two  accelera¬ 
tions  will  act  on  the  point  mass  missile.  The  sum  of  the  linear  accelera¬ 


tions 


_ 

R  =  g  +  a. 


=  g  +  a,p| 


(6.1) 


where 


R  =  resultant  inertial  acceleration  of  the  missile  C..G. 


a,p  =  magnitude  of  thrust  aoceleration  (=  thrust  per  unit 
mass) 

I  =  unit  vector  in  the  direction  of  the  thrust  vector 

^ 

The  ballistic  missile  must  be  steered  such  that  R  -'■Vt,  (R,  t)  at  burn- 

X\ 

out  for  nominal  as  well  as  nonstandard  missiles.  The  velocity  to  be 

gained,  V  =  V  (R,  R,  t)  is  defined  as 
g  g 


=  Vj^(R,t)  -  R 


(6.2) 


Hence,  a  commanded  thrust  attitude  vector,  (t),  must  be  found  such 
that  all  three  components  of  vanish  at  burnout.  At  the  precise-  time 
when  Vg  =  0  the  missile  engine  will  be  cut  off.  The  missile  will  then 
coast  to  its  target  satisfying  the  four  guidance  restraints. 

It  is  convenient  to  choose  the  computing  coordinate  system  set  of 
basis  vectors  as  follows: 

i  =  parallel  to  the  nominal  thrust  attitude  vector  at 
burnout 
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j  =  Pj^,  i.e.,  normal  to  the  pitch  plane 

"k  =  TaT 


The  origin  of  the  inertial  coordinate  system  may  be  chosen  at  the  center 
of  the  earth  or  at  the  launch  site.  The  velocity  of  the  missile  will  be 
^ 

R  =  Xi+Yj+Zk 

s 

Pitch  commands  may  then  be  based  on  Z^,  yaw  commands  may  be  based 

on  Yg,  and  engine  cutoff  commat^s  based  on  X^.  These  components 

are  defined  by  equation  (6.2),  and  Vj^  is  obtained  by  explicit  computation 

ot  by  using  one  of  the  methods  given  in  Section  V.  It  will  be  assumed 

that  the  change  in  thrust  attitude  is  sufficiently  small  that  cross -coupling 

between  X  ,  Y  ,  and  Z  may  be  neglected.  „ 
g  S  g 

When  the  thrust  attitude,  is  specified,  the  motion  of  the  rhissile 
C.G.  is  obtained  by  integrating  equation  (6.  1).  Since  4  is  a  unit  vector 
the  specification  could  be  given  in  terms  of  the  pitch  and  yaw  components, 

£  and  £  .  The  specification  could  also  be  given  in  terms  of  the  pitch 

*  y  ♦  • 

and  yaw  attitude  rates,  £  and  £  .  It  is  assumed  that  the  commanded 

z  y 

roll  attitude  rate  is  zero.  The  following  paragraphs  will  discuss  the 
determination  of  the  pitch  attitude.  A  corresponding  discussion  (not 
presented)  will  be  obvious  for  the  yaw  channel. 


The  pitch  attitude  of  the  thrust  vector,  £  ,  must  be  determined 

z 


such  that  V 


0.  This  result  would  be  obtained  if  V  were  a  solution 


of  any  one  of  the  following  six  forms; 


1)  -Z  +  KV 


V  +  KV 
gz  gz 


3)*  -Z  +  KV 


It  is  understood  throughout  that 
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0 


(6.4) 


4)  V  +  KV  +  L  r  V 

gz  gz  j  g® 

to 

5)  -Z*  +  KV  +  LV  =0 

gz  g* 


6)  V  +  KV  +  LV  =0 
gz  gz  gz 


•where  =  ^Rz  ~  other  forms,  including  cross -coupled  forms, 

could  be  added  to  the  above  list,  but  the  above  six  will  suffice  for  the  pre¬ 
sent  discussion.  The  gain  functions  K  and  L  are  chosen  to  insure 
stability  and  satisfactory  transient  response.  The  presence  of  the  V^^ 

term  in  all  six  forms  assures  that  V  will  be  controlled  to  zero. 

gz 

Forms  (1)  and  (2)  are  the  simplest,  form  (2)  being  preferred  because 
Vg^  will  also  vanish  at  burnout.  In  practice,  it  is  found  that  the  presence 
of  thrust  misalignments,  lags  in  the  control  system,  etc. ,  result  in 
steering  errors.  These  errors  may  be  significantly  reduced  by  intro¬ 
ducing  integral  control  terms  as  in  forms  (3)  and  (4).  In  these  forms, 
t^  designates  the  time  at  which  guidance  is  initiated.  Forms  (5)  and  (6) 
are  equivalent  to  forms  (3)  and  (4)  and  are  used  when  the  thrust  attitude 
rate  is  to  be  commanded. 

•  • 

When  form  (4)  is  used,  ^  is  obtained  by  eliminating  Z  using 
equatipn  (6.  1)  as 


Z  = 


gz  *  ^T^z 


(6.5) 


Then 

a_|  =  KV  +  L  V  +Vo  -g  (6.6) 

T  z  gz  gz  Rz  «z 

Equation  (6.6)  may  be  simplified  when  V„  (R,  t)  is  chosen  from  a 
restricted  class  of  vector  functions.  The  restriction  r^uires  that  when¬ 
ever  the  missile  engine  is  shut  off,  with  V  =0,  that  V  will  continue 

g  S 

to  vanish  during  the  coast  period  that  precedes  arrival  at  the  target. 
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Guidance  constraints  of  the  type  a)  -  c)  on  Page  21  will  not  result 
in  expressions  for  Vj^  (R.  t)  that  satisfy  the  class  restriction.  Con¬ 
straints  d)  -  h)  will  satisfy  the  restriction.  In  general,  the  class  restric¬ 
tion  is  not  a  serious  limitation.  For  example,  whenever  the  four  guidance 
restraints  are  defined  at  the  target,  the  restriction  will  be  satisfied.  The 
remaining  analysis  in  this  report  will  assume  that  Vj^  satisfies  the  class 
restriction. 

-A  ^ 

When  the  guidance  restraints  are  chosen  such  that  Vj^  (R,  t)  satisfies 
the  class  restriction,  it  follows  that 


g  =  Vj^  (R.t).  Vp^  =  R 


(6.7) 


The  following  expansion  in  partial  derivatives  is  valid  in  general 


a=l 


(6.8) 


where 


Rj  =  X,  R2  =  Y,  and  =  Z 


When  Vp  satisfies  the  class  restriction,  equations  (6.7)  maybe  substi¬ 
tuted  into  (6.  8)  to  obtain 


=  y  V  +-^ 

Ra  St 


(6.  8A) 


where 


^R1  ""  ^Rx’  ^R2  “  ^Ry’  ^R3  “  ^Rz 


Elimination  of  — from  equations  (6.8A)  and  (6.  8)  gives 


^R  "Z- 


(6.9) 
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This  simplification  will  not  be  made  in  this  section  because  of  the  com¬ 
plex  character  of  equation  (6.  11).  However,  it  is  felt  that  constants  K 
and  L  for  the  pitch  (and  yaw)  channels  may  be  fo\md  that  will  cause 
V  -*  0  and  V  -*■0.  Then  the  V  terms  in  equation  (6.  10)  will  not 

gz  gy  gQ  H  \  I 

be  required. 

When  a  missile  is  steered  according  to  equation  (6»  10)  an  excessive 
use  of  propellant  may  be  experienced,  since  the  nominal  trajectory  will 
not  be  reproduced.  The  nominal  pitch  thrust  attitude  is  obtained  as 

^Nz  "  ^ 

If  ^  were  evaluated  along  the  nominal  trajectory,  it  would  generally  be 
z 

found  that 

^z  ^  ^Nz 

By  construction  of  it  follows  that  any  such  deviations  can  only  result 
in  the  failure  of  equation  (6.  10)  to  satisfy  the  nominal  trajectory  shaping 

—A 

restraints,  with  the  exception  of  the  first  (R-*  Vj^).  This  defect  in  equa¬ 
tion  (6.  10)  can  be  remedied  in  two  ways,  both  essentially  equivalent.  A 
discussion  of  this  procedure  follows. 
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The  reason  why  equation  (6.  10)  causes  the  nominal  missile  to  be 
steered  (in  pitch)  is  that  equations  (6.4)  are  not  satisfied  by  the  nominal 
missile.  Equations  (6.4),  however,  can  be  easily  generalized  so  that 
they  are  satisfied  by  the  nominal  missile.  This  is  done  by  adding  a 
forcing  function,  F  (V  ),  to  the  right-hand  sides.  The  forcing  function 
is  obtained  by  evaluating  the  left-hand  side  of  equation  (6.  4),  form  (4), 
over  the  nominal  trajectory,  to  obtain  as 


V  . .  +  K  V  . . 
gzN  gzN 


(6.  12) 


gxN’ 
The 


The  nominal  value  of  the  X-component  of  the  velocity  to  be  gained,  V 

N 

is  also  obtained.  Then  F^  is  expressed  as  a  function  of 

quantity  is  chosen  as  the  independent  variable  for  the  F-function 

because  V  will  tend  to  zero  in  a  smooth  fashion  as  the  missile  gains 
gx  ® 

=  0;  hence,  0  as  required.  In  order  to 


velocity.  At  burnout  V 
retain  computer  simplicity,  F^^  is  chosen  in  the  form 


gx 


N 


gz 


.N 


=  C,  V  „  +  C,  V 

1  gxN  2  gxN 


(6.13) 


where  and  are  trajectory  shaping  constants.  These  constants 

should.be  obtained  such  that  a  good  fit  to  the  nominal  trajectory  is  obtained 
in  the  burnout  region.  The  curve  fit  should  also  include  the  guidance 
initiation  point  in  order  to  minimize  pitch  commands  when  the  guidance 
loop  is  closed. 

Although  quadratic  terms  are  indicated  in  equation  (6.  13),  in  some 
cases  linear  terms  will  suffice;  in  others  cubic  terms  will  be  called  for. 
This  choice  depends  on  the  precision  required  in  duplicating  the  nominal 
trajectory.  Only  a  small  amount  of  propellant  waste  is  allowed  for 
missions  requiring  maximum  burnout  energy.  In  such  missions  every 
additional  X  pounds  of  propellant  required  represents  the  loss  of  an 
additional  physical  experiment. 


The  forcing  function,  added  to  the  right-hand  side  of 

equation  (6.4),  fornm  (4)  is  then 


F  (V  ,)  =  C,  V  ^  +  C,  V 
z  '  gz  1  gx  2  gx 


(6.14) 


where  Cj^  and  are  the  same  constants  used  in  equation  (6.  13). 

Equation  (6.4),  form  (4),  is  then  modified  as 


V  +  K  V  +  L 


'  =  F  (V  ) 

gz  z  gx' 


(6. 15) 


Equation  (6.  15)  is  now  satisfied  by  the  nominal  missile, and  V  -►  0  as 

a 

required.  The  yaw  and  pitch  attitudes  to  be  commanded,  ^  and  |  , 

y  z 

may  now  be  obtained  following  the  procedure  indicated  by  equations  (6.  5) 
to  (6.  10).  The  result  is 


di-i  =  K  V  +  L 

T*y  gy 


a_e  =  K  V  +  L 
T^z  gz 


L 

r‘v 


SE  V  -  F  (V  ) 
Q  y  '  gx' 


(6. 15A) 


V  -  F  (V  ) 
9X0  ga  z  '  gx 


where  it  is  assumed  that  the  pitch  and  yaw  channels  are  identical;  hence, 
the  same  gains  K  and  L. 

A  second  procedure  will  now  be  mentioned  for  shaping  equations 
(6.4)  so  that  these  equations  are  satisfied  by  the  nominal  missile.  Again 
let  equation  (6.4),  form  (4),  be  used  as  the  example. 


Note  that  equations  (6.6),  (6.  10),  (6.15A),  and  (6, 18)  can  be  simplified 
if  the  computing  coordinate  system' is  chosen  such  that 


Vp  -  g  /si  =  0 
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The  desired  forcing  function  is 


'  '"gzN  +  K  ’^gzN  +  ^  /  '"gzN 

‘^to 


Hence,  form  (4)  could  be  written  as 


•  #  «  r'-  * 

V  +KV  +L  V  =  0 
gz  gz  /  gz 

“'to 


(6. 16) 


where 


V  =  V  -  V 
gz  gz  gzN 


(6.17) 


The  quantity  may  now  be  expressed  in  powers  of  as  was 

done  in  equation  (6.  13)  for  F^.  Note  that  this  method  requires  that 
fitting  errors  in  the  integral  and  the  derivative  of  also  be  mini¬ 

mized,  When  this  technique  is  used  the  expressions  for  and 
are  obtained  as 


a_5  =i_6  *  KV  *  +  L  f  V* 

TV  T_^V  gy  gy 


0V 

y  V* 

,  (6.18) 


a„^  =  a„  I  +  K  V  +  L 

T’z  T^  zn  gz 


f  V*  -  y 
J  to  9 


V* 

ga 


^  ^  9Q(  sje 

where'  V  ,  =  V  and  V  _  =  V  are  defined  as  in  equation  (6,  17). 

gl  gx  g2  gy  'i  ^  f 

That  is 


V  =  V  -  V  . , 
g  g  gN 


(6.  19) 


Since  V  -*•  0  at  burnout,  it  follows  that  the  variation  between  a,,,  and 
g  T 

a,pj^  may  be  neglected.  Then  equations  (6. 18)  become 
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=  e  JL  V*  +  i  r‘  V* .  i.  r 

yN  gy  gy  2^ 

+  Jl  f\*  --1  y 

a*  /  gz  a  Z. 
-'to 


0V 


^  V* 


8X‘‘  8® 


(6.20) 


£  +  V* 

^zN  ^  a^  ''gz 


8V 


Rz  y* 


0X“  80 


Equations  (6.20)  completely  determine  since  it  is  only  necessary  to 
specify  two  components  of  a  unit  vector,  because  equations  (6.  20)  could 
also  be  written  as 


4^ 


e 


K  * 

+  #-v 

gy 


K  * 

+  — V 

gy 


(6.20A) 


where  ij/  -  and  0  -  0j^  denote  small  variations  in  the  pitch  and  yaw 
attitude  angles  about  the  nominal  values. 


Equations  (6. 15 A)  and  (6.20)  may  now  be  used  to  steer  the  ballistic 

missile.  The  quantity  will  furnish  the  engine  shutoff  command. 

That  is,  when  V  =  «  a  command  will  be  issued  to  shut  off  the  main 

gx 

engine.  The  parameter  €  is  included  to  allow  for  the  expected  residual 
impulse  of  the  main  engine,  as  well  as  the  impulse  of  a  vernier  engine, 
if  provided.  The  vernier  engine  will  be  shut  off  when  V  =0. 

9it  ^ 

Since  V  ^0,  equations  (6.  20)  can  be  written  in  vector  form  as 

gx 


’N  a„  g  a„ 


.±  y 

a™  L. 


9V 


to 


0X 


a  ga 


(6.21) 


where  the  gain  constants  K  and  L  are  also  indicated  for  the  roll 
channel.  This  is  done  simply  for  mathematical  convenience  and  is  accept¬ 
able  since  it  is  not  planned  to  steer  in  roll.  The  symbol  is  used  in 
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equation  (6.21)  to  denote  parallelism.  The  X  axis  was  chosen  along 


at  burnout.  Hence,  if  is  a  constant  vector,  it  follows  that  equa¬ 
tion  (6.  21)  is  normalized  except  for  second-order  terms.  When  is 
not  a  constant  vector,  it  is  desirable  to  define  the  quantities  such 
that 


it 


=  V 


=  A 


•N 


g 


’N 


(6.  22) 


By  construction,  vt  is  the  sum  of  the  two  components  of  V  normal  to 

g  ^  g 


the  thrust  attitude  vector  ^ 


N* 


will  be  denoted  by  e,  hence 


The  component  of  V  parallel  to 

g 


Then 


g 


^g 

=  6  i 


'N 


TVT  +  vt 

N  g 


(6.23) 


Both  V  and  V'  vanish  along  the  nominal  trajectory  and  are  identical 
if  is  a  constant.  Hence,  equations  (6.  21)  may  be  written  as 


f  =  -S.  vt  4  i  f‘vi.±y 

N  a„  g  a„  I  g 

^  ^  to  ^  1 


3  8V 


8X 


5:  ^t 

a 


g“ 


(6.24) 


Equation  (6.  24)  may  also  b^  used  for  missile  steering.  The  main 
engine  cutoff  signal  is  e,  defined  by  equation  (6.  23).  The  advantage  of 
this  representation  is  that  it  is  independent  of  a  particular  coordinate 
system  and  is  free  of  cross -coupling  errors. 

The  thrust  attitude  rate,  |,  may  be  commanded  in  place  of  This 
expression  may  be  derived  by  writing  equation  (6.4),  form  (6),  as 


Vg  +  K  Vg  +  L  Vg  =  q  (e) 


(6.25) 
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The  inhomogeneous  term  q  (e)  is  added  so  that  equation  (6.25)  will  be 
satisfied  by  the  nominal  trajectory.  That  is 


(6.26) 


where  e  denotes  the  engine  cutoff  signal.  Equations  (6.  1)  and  (6.  2) 
may  be  differentiated  with  the  result 


R  =  g  +  a.j.  ^  +  a.j.  ^ 


Using  the  above  results,  equation  (6.25)  becomes 


’T4a«  4  = 


Vo  -g  +K  V„  +  LV„  .q(e) 


(6., 27) 


The  quantity  a.^.  g  may  be  absorbed  in  q  (e).  Then  equation  (6.  27)  will 

specify  a  commanded  turning  rate  that  will  c;au8e  V  -*0.  Since, 

g 

V  =  V_  -t.a,^t 


eqviation  (6.  27)  may  be  written  as 


aJ*  +  +  Ka„)  t  = 


L  V  +  K  (Vj^  -  g)  +  (Vj^  -  g)  -  q  (e)  (6.28) 


The  result  of  multiplication  by  ^  (vector  cross  product)  is 


,\Vg  +  e,^  Vj^-g  + 


K  ^  A  Vp  -  g  -  6  A  q  (e) 


(6.29) 


where  u  is  the  commanded  angular  velocity  of  the  thrust  attitude  vector 
defined  as 


Q  =  I  A  I 


(6.30) 
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1 


The  steering  equation  (6.29)  is  one  form  of  the  so-called  cross-product 

steering.  This  equation  may  be  simplified  using  equations  (6.  9)  and 

choosing  the  coordinate  system  suggested  in  the  footnote  on  Page  37. 

Other  simplifications  will  result  when  suitable  cross -coupled  terms  are 

added  to  equation  (6.4),  form  (6)  as  discussed  in  connection  with  equation 

(6.  11).  A  stability  analysis  should  then  be  done  on  the  resulting  form. 

Finally  the  steering  equations  should  be  simulated. using  a  closed-loop 

guidance  simulation.  This  last  step  is  necessary  to  verify  that  only  . 

negligible  errors  (V  ^0  'at  burnout)  result  when  nonstandard  missiles 

S 

are  steered. 
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VII.  BALLISTIC  MISSILE  CONTROL  SYSTEMS 


The  preceding  discussion  was  restricted  to  point  mass  missiles, 
where  the  thrust  chamber  was  located  at  the  C.G.  Since  this  is  not  the 
case  for  long  slender  ballistic  missiles,  it  is  necessary  to  consider 
actual  missile  dynamics.  Only  rigid-body  missile  dynamics  will  be  con¬ 
sidered  in  this  section.  Aerodynamic  forces  will  nbt  be  treated  since  it 
is  assumed  that  the  missile  is  above  the  atmosphere.  The  forces  acting 
on  the  missile  are  shown  in  Figure  7.1.  In  Figure  7.  1  the  thrust  attitude 
unit  vector  is  denoted  by  the  symbol  6.  In  the  preceding  section  the 
thrust  attitude  was  denoted  by  This  symbol  will  now  be  used  to 
denote  the  missile  attitude  vinit  vector.  When  this  change  in  notation  is 
made,  the  linear  acceleration  equation  becomes: 

^  =  g+a,j.6  (7.1) 


Two  additional  unit  vectors  are  defined  by  the  engine  missile  gimbal 
geometry.  Pitch  rotations  of  the  missile  occur  about  yaw  rotations 
about  "q. 

The  relation  between  missile  attitude,  and  thrust  attitude  will 
now  be  found.  This  determination  requires  application  of  the  turning 
moment  equation.  The  moment  equation  for  missiles  is 

Iw  (7.2) 

torque  about  the  missile  C.  G. 

I  (t)  is  the  missile  pitch  or  yaw  moment 
of  inertia  about  the  C.  G. 
d 

ar 

e  At 

Equation  (7.2)  neglects  propellant  sloshing  effects;  it  is  assumed  that  all 
internal  linear  momenta  are  directed  along  the  missile  axis.  In  order  for 
equation  (7.2)  to  have  the  given  simple  form,  the  roll  angular  velocity 
must  vanish.  Then  the  roll  moment  of  inertia  is  not  involved.  The  torque 
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where 


a  = 

I  = 

«  = 


w  = 


Figure  7.1  -  Forces  Acting  on  the  Ballistic  Missile  During 
Powered  Flight  in  a  Vacuum. 


about  the  missile  C.G.  is 

5“=  T  ^  (it) 


i  t6 


Hence , 


u  6  A  'fe  =  ci 
c 


where 


(7.3) 


(7.4) 


(7.5) 


Equation  (7.4)  is  usually  presented  in  component  form.  The  pitch  compo¬ 
nent  is  obtained  by  multiplying  (scalar  product)  equation  (7.4)  by  (  to 
obtain 

p  6  %'V 

'^c  o 


where 


6^  =  pitch  angle  between  6  and  ^ 


0 


%  6  ^  •  C 

magnitude  of  pitch  angular  acceleration 


~  u> 


The  symbol  %  is  used  above  to  indicate  that  small  angle  approximations 
have  been  made. 

The  purpose  of  the  missile  control  system  is  to  steer  the  thrust 
vector,  T,  such  that  the  missile  will  assume  a  reference  attitude,  | 

The  reference  attitude  is  chosen  to  coincide  with  the  desired  thrust  attitude 
determined  in  the  preceding  section.  The  missile  moment  equation  (7.4) 
indicates  the  validity  of  this  identification.  When  ^  ^  6,  the  missile  must 
experience  an  angular  acceleration,  w.  The  assumption  that  u  is  negli¬ 
gible  is  acceptable  when  guidance  steering  commands  only  are  considered. 
That  is,  u  is  large  as  a  result  of  control  system  commands  rather  than 
guidance  system  commands.  When  w  =  0,  equation  (7.1)  reduces  to 


R  =  g  +  I 
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This  is  precisely  the  same  equation  that  was  discussed  in  Section  VI. 

Expressions  for  the  reference  attitude  were  given  in  equations  (6.  1^), 
(6.18).  (6.20).  (6.21).and  (6.24).  The  latter  is 


t  +  f* j.y 

N  a.p  g  a.j.  J  g  a_^ 


The  commanded  thrust  attitude.  6  ,  must  now  be  determined  to  insure 
that  the  missile  attitude  will  be  controlled  to  . 


The  commanded  pitch  or  yaw  thrust  attitude  is  usually  presented  in 
component  form,  the  pitch  component  (typical)  is 

«oc  =  (Sr  -  «>  - «  (■'•■'> 

where 

6  =  pitch  component 'of  commanded 

oc  \  s, 

thrust  attitude  6^  •  O 


Kj^.  Kj^  =  control  system  constants 

—X  — i  \ 

-  0  =  pitch  angle  between  and  C  (  %  ^  A 

'  — >  — i 

0  =  magnitude  of  the  pitch  angular  velocity  (%a)  .  C) 


In  equation  (7.7)  the  quantity  0  is  usually  measured  by  a  rate  gyro  pack¬ 
age  located  at  some  point  in  the  missile  structure  that  is  relatively  free 
of  bending  frequencies. 

The  control  system  equations  may  be  specified  in  a  form  that  is  in¬ 
dependent  of  any  particular  coordinate  system.  Equation  (7.  7)  is  readily 
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(7.8) 


generalized  to  the  vector  form.  The  result  is 

Equation  (7.8)  defines  both  pitch  and  yaw  engine  deflection  commands. 

The  same  constants  Kj^  and  Kj^  are  used  for  both  pitch  and  yaw  since 

these  channels  are  assumed  to  be  identical.  In  this  section,  as  in  the 

preceding  section,  the  lag  in  the  engine  hydraulic  system  will  be  neglected. 

Hence,  the  quantities  6  =  '6^  will  be  considered  to  be  identical. 

When  6=6  is  eliminated  from  equations  (7.8)  and  (7.4)  the  result  is 
c 


ui  +  A  %  0 


(7.9) 


The  pitch  component  of  equation  (7.9)  is  obtained  by  multiplying  (scalar 
product)  equation  (7.9)  by  The  result  is 


ot^^Kj,o  +  ^^Kp(e-ej,)so 


(7.  10) 


Equation  (7.  10)  is  usually  obtained  by  eliminating  6^  =  6^^  between 
equations  (7.  7)  and  (7.6).  A  consideration  of  equation  (7.  10)  serves  to 
justify  the  form  chosen  for  6^  in  equations  (7.  7)  or  (7.8).  The  control 

system  constants  and  K  must  be  chosen  such  that  a  stable  system 

L)  jK 

with  satisfactory  transient  response  is  obtained.  In  practice,  missile 
bending. and  propellant  sloshing  must  be  considered  so  that  a  simulation  is 
required  for  this  gain  determination.  Values  of  and  Kj^  may  gener- 

ally  be  found  such  that  t  **  with  small  lag. 

If  some  additional  computer  complexity  is  allowed,  then  can  be 
calculated  as  in  equation  (6.27).  The  pitch  component  (typical)  of  6^ 
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could  then  be  written  in  the  form: 


Equations  of  this  type  wotild  control  5  to  as  well  as  t  to  6^,  a 
desirable  feature  when  the  nominal  missile  is  pitched  over  at  a  high  rate. 
Some  additional  relations,  useful  in  guidance  and  control  systems 

-i. 

analysis,  will  now  be  given.  Since  ^  A  f  ,  equation  (7.  8)  may  be 
written  as: 


6c  A  e  =  Kjjg  -  Kj^  eAe 


(7.11) 


Equation  (7.  11)  will  now  be  assumed  to  define  6c,  hence.  the  omission  of 
the  symbol.  An  equivalent  form  of  equation  (7.  11)  is 


I  A 


«c  V 


=  0 


(7.12) 


The  quantity  in  the  square  parenthesis  must  be  parallel  to  i,  hence, 

Ae  =  6c  +  -  Kj^l  (7.  13) 

The  normalization  constant.  A,  may  be  determined  by  multiplying  (scalar 
product)  both  sides  of  equation  (7.  9)  by  |  .  Then 


A  =  6c  .  e  +  .  i 


(7. 14) 


The  result  of  substituting  A  into  equation  (7.  13)  is 

t  =  t  (t  .  t)  -  K 


D 


+  K^t  (7.15) 


When  the  small  angle  approximation  is  made,  equation  (7.  15)  becomes 


(7. 16) 
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^  —A 

Equation  (7.  16)  gives  a  simple  relation  between  6  ,  and  ,  that 

C  n. 

does  not  depend  on  any  particular  inertial  coordinate  system. 

— i  —Ik 

A  vector  relation  may  also  be  found  relating  6  and  &  .  Eqtiation 
(7.4)  is 


-A 

=  d  (C  a4) 
It 


(7.17) 


Since  =  0  ,  the  moment  equation  reduces  to 


“c'  *  i 


(7.18) 


It  may  be  shown  that 


-  + 


w 


1 

b> 


(7.  19) 


where  «  and  u  are  the  magnitudes  of  w  and  ~w.  With  this  substitution, 
equation  (7.  18)  becomes 


eA 


6  + 


b) 


u 


=  0 


(7.20) 


An  expression  for  6  may  be  obtained  from  equation  (7.  20)  using  the 
method  of  the  preceding  paragraph.  The  result  is 


6  = 


e)-  ^ 


(7.21) 
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Since  i  and  ^  are  unit  vectors,  it  follows  that 

h> 

cos  6j  =  6  .  6 

.  w 
sin  6,  =  — 

1 


(7.22) 


where  6^  is  the  angle  (combined  pitch  and  yaw)  between  if  and  When 
the  small  angle  approximation  is  made,  equation  (7.  21)  becomes 


^  fc  ~4  -  —  i 


(7.23) 


The  above  discussion  presented  the  basic  equations  of  a  missile 
control  system  necessary  for  guidance  system  analysis.  The  discussion 
assumed  that  the  missile  attitude  was  to  be  controlled  to  the  guidance 
reference  attitude  .  Only  simple  modifications  are  required  to  adapt 

^  -1..  -J.  -A 

the  analysis  to  the  case  when  |  is  to  be  controlled  to  or  u  to  . 

A  small  complication  could  be  added  to  the  above  analysis  by  includ¬ 
ing  the  effect  of  the  engine  gimbal  inertia.  The  relationship  between  6 

—A  k 

and  6^  is  used  in  this  section  as  6  =  6^  .  A  better  approximation  for 
this  relation  is 


*o<‘> 


-kr 

6  (t  -  t)  dr 
oc 


(7.24) 


for  the  pitch  (typical)  component.  In  equation  (7.24)  the  parameter  k  is 
related  to  the  lag  in  the  engine  hydraulic  system. 
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APPENDIX  A 


Derivation  of  the  Force  Equation  for  Ballistic  MiBsileg 

A  short  nonrigorous  derivation  of  the  linear  force  equation  is  given 
in  this  appendix.  The  sum  of  the  .external  forces  acting  on  a  system  of 
particles  equals  the  rate  of  change  of  linear  momentum  of  the  system  or 

F  =  (mv  +  m  V  )  (A-1) 

^  g  g 

where 


m  =  mass  of  the  missile 
V  =  velocity  of  the  missile 


m  =  mass  of  the  escaping  gas 
g 


V  =  velocity  of  the  escaping  gas 
g 


and  the  dot  is  used  to  denote  differentiation  with  respect  to  time.  It  is 
assumed  in  this  appendix  that  the  only  external  force  on  the  missile  arises 
from  gravitational  acceleration;  hence,  equation  (A-1)  may  be  written  as 


mg  =  mv  +  mv+mv  + 


m  V 
a  cr 


(A-2) 


Now  V  =  0  when  the  gas  exits  into  free  space,  and  m 
8 

total  system  mass  is  a  constant.  Thus 


=  -m  since  the 
g 


where 


C  =  escape  gas  velocity  with  respect  to  the  missile 
(also  called  specific  impulse) 

T  =  -  m  C  =  missile  thrust  vector, 

g 

Equation  (A- 3)  may  be  divided  by  m,  with  the  result 

R-  =  "v  =  "g  +  "a^  (B-4) 

where 

-X  T 

~  in" 
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APPENDIX  B 


The  Gravity  Tarn  Pitch  Program 

The  missile  is  steered  through  the  atmosphere  such  that  a  "gravity 
turn"  is  followed.  This  pitch  profile  is  alternately  referred  to  as  either 
a  zero -lift  or  a  zero  angle -of -attack  pitch  program.  This  program  is 
utilized  in  order  to  prevent  breakage  of  the  missile  as  a  result  of  aero¬ 
dynamic  forces. 

The  aerodynamic  forces  are  referred  to  as  drag,  D,  and  lift,  L, 

The  drag  force  is  directed  along  the  roll  axis  while  the  lift  force  acts 
normal  to  the  missile  axis.  The  forces  acting  on  the  missile,  for  flight 
through  the  atmosphere,  are  shown  in  Figure  B  -  1^,  The  point  of  applica¬ 
tion  of  the  resultant  aerodynamic  force  is  referred  to  as  the  center  of 
pressure  (C.P. ).  The  forces,  L  and  D,  acting  at  the  C,  P,  are  defined 
as; 


Figure  B-1  -  Forces  Acting  on  the  Missile  for 
Flight  Through  the  Atmosphere 
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^  .A 

D  =  -  D  4 
L  =  L  tA(V.At) 

where: 

Cd  A  q 
Cj^  ia)  A  q 

direction  of  missile  roll  axis,  a  unit  vector 

velocity  of  the  missile  \vith  respect 
to  the  air  mass 

magnitude  of 

effective  missile  area 
2 

1/2  p  =  dynamic  pressure 

air  density 
drag  coefficient 

lift  coefficient 

angle  of  attack  cos  a  =  ■  £) 

The  drag  and  lift  forces  are  sometimes  defined  as  acting  along  and 
normal  to  V^,  rather  than  |  When  actual  data  are  used,  care  must  be 
taken  to  insure  that  the  definitions  are  consistent  with  the  data. 

The  axial  strength  of  the  missile  is  greater  than  the  transverse 
strength.  Hence,  the  normal  forces  (lift)  must  be  minimized  for  flight 
through  the  atmosphere.  Otherwise  the  aerodynamic  lift  forces  would 
produce  bending  moments  that  could  break  the  long,  slender  missile. 


D  = 
L  = 

t  = 


A  “ 
A  = 

q  = 

p  = 
'D  “ 


CN(a)  = 


(B-1) 

(B-2) 
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Note  also  that  when  L  ^  0,  it  is  necess-ary  to  choose  the  engine  thrust 
direction,  6,  to  prevent  the  missile  from  rotating.  The  aerodynamic 
pitching  moment  is  cancelled  by  choosing 

Ld  =  fA{6At)^  T 


=  [t-tit  .  6)]  Tl 


where 


d  =  distance  between  the  C.G.  and  the  C.P. 
f  =  distance  between  the  C.G.  and  the  engine  gimbal 
T  =  engine  thrust 


When  a  zero -lift  pitch  program  is  not  followed  the  energy  required  to 
cancel  the  pitching  moment  is  wasted.  Some  of  this  wasted  energy  is 
converted  to  heat  energy  resulting  in  weakening  of  the  missile  structure. 
Note  that  L  can  be  made  to  vanish  by  choosing 


i  = 


(B-3) 


The  quantity  is  calculated  as 


■  B.-Wj,AR+W 


(B-4) 


where 

R  =  velocity  of  the  missile  with  respect  to 
inertial  space 

R  =  position  vector  of  the  missile  with  respect  to  the  earth 

=  angular  velocity  of  the  earth 

W  =  wind  velocity  with  respect  to  the  earth, 
normally  a  negligible  quantity 
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The  missile  will  fly  a  gravity  turn  when  equation  (B-3)  is  satisfied. 

-A 

All  the  quantities  required  to  determine  the  direction  of  are  calculated 

in  standard  missile  simulations.  Then  the  thrust  attitude,  6,  may  be 
commanded,  as 


When  aunity  control  system  (equivalent  to  a  point  mass  missile)  is 
assumed,  it  follows  that  6  =  ^. 


The  simulation  will  calculate  the  missile  angular,  velocity,, 
cj  =  ^  A  This  quantity  may  be  approximated  with  simple  ftmctions 
and  incorporated  in  a  missile  pitch  programmer.  Steering  commands 
from  the  missile  programmer  then  cause  the  missile  to  pitch  over  to 
approximate  the  desired  gravity  turn. 

Note  that  when  W  %  0,  the  initial  value  of  |  =  6  is  arbitrary. 
’ATien  the  initial  value,  ^  is  chosen  the  resulting  gravity  turn  is 

specified.  In  general  it  is  necessary  to  make  several  flights,  using 
different  values  for  4^,  in  order  to  obtain  the  desired  end  conditions. 

These  variations  can  be  made  even  when  W  ^  0,  since  wind  velocities 
are  normally  small  enough  to  be  neglected. 


-57- 


2  0 

.a  R  J  -  7  * 

5 wa- K  -0 
a 

'S 

S 

|3<2S;  I 

Isiig 

“illsj 

IgHSSK 

0  e  j  ^  e<  5 

frga  .-e 

i>  S  Bd  'A  ^  f^j 

a§o^ 

•  g'^b  a  2 
o  Hwgj  a  • 

Jj  2  H  O  -  6 
V w Q, a,  .a 


^5s3i 


•d  **  S 

h  II  >;  II 

**  •**  m  ^ 

5  S-jg^  ?« 

g5  l-g 
a  s  1  s 

5  C  II  'j  s 

^  •S'i  B  3 

g  B  2  Jl 

u  •  •!  i«  6 

B  ii  J  fcl 

'*  "S  j  w  ^ 

«  ^  5  I) 

S? 

I  •»  S  g  Jj  >0 

“  B  ■-  C  II  u 

M  11  “jj  a  g 

6  •  •  5  —  fi 

•  £  *  M  '3  8 

liflli 

*•2*25  S 

JV  S  £  0*  M  h 

>  E  ?  9.3  I* 


j3  0  ■  2 

g  °  JS  .£  - 

"  S’g'S  ?•« 

•2|^  So-S 

T>  7,  o  0.2  -s 

J  U  *  0 
U  g  X  2  i!  II  a 

§:«  S  gi-S- 

**  I  ir'?'’S  S 

a  .«  a  ea-g  « 

a  *  sa  a 
•  aj  1  "If 

•3  0  5  5  o  § 

u  *  •  B  'S 

« '21  is  I 

wjr  II fH  •  ^ 

22II,.| 

S  a  •?  II  ’c.  a 
m~  a  M  “  0 
5  S'  2  B  *>  U 


2  o 

B  R  J  "  r*  5 

B  ^  .i  a 

g«2-H  "S 
a  5  H  o  S  “ 

B  <  g  a  g:  B 

^Sh5-  5 
|g<2S?  I 
I^S  bQ 

Ji  H  >3  JJ  H 

71“*”^  • 
w  W  g  O  o 

.22  ii2X 

gw  w  ^ 

o  o  ^  *r  H  ^ 

a  2  n  ^  0  “? 

2 »p 

o  ^  <  ci  g  B 

5gg-»3^ 

S-SSp-gt;- 

o  H  v)  s«  C  2 
2  H  O  ^ 

mu 


w  a  ” 
u  •>  TS  JJ  £ 

.H  M  ^  V 
^  ‘*B  ■  _P 

;3  S-^B  J. 

S  t)  •"  o  « 
l«  M  •  O  M 

•“  B-S  £.£ 

m3  a  “  * 

£  ■“  E  .5  E 
-■«  §  5.S 
£  B  SI'S 
b  -S’3  b  * 
o  2  2  B  « 

o  5  B  a  6 

^  SS  2l 

4)  j  ^  *1 

O  <  .H  .  y  a 

•3  «•  B  Ji  £  “ 
£?  t-S 

&B-|E  «  B 

w  4)  -  >4  e 

®  *  W  "S  w  ^ 

"  Q.  K  B  a  • 

u  “■  O  “  M  6 

V  ft  «4  ^  I)  W 

g*  44  ij  c  •g 

5:5-«£ 

II  si?; 

ill  15  s 


4!  O  *0  •  £ 
g  w  4>  V  h 

3  w  .bb  ^ 

•  J  O  ^  ?-tJ 
B  M  a  B  .5  ^ 

■"  3  B  £  u  3 

■o  7.  o  0.2  -a 

t>  g  o  *0  9 

S  g  •  o  oe^ 
U  g  S  2  £  B  a 

g(J:i  °b|s 
S'  .1  B^-s-^ 

o.  fi  '  o  ■“  5 

B  i i>>.7  £ 

B  J3  a  *  B  u 

a  -S  a  0''S  *> 

B  P  S  «  .S  * 

«  a.B  Saf 

S'"  'll  o  s 

■fi  I  B  *  O  -2 

o  ,  .2  £•  B  u 

2  -  E  0 .2  a 

.Sa«  S2  ? 

a  3  - 1 H  -a 

£33  •  I 

S  -  ®  .2  B  •£ 

•  S  B  O'  >  S 
!h  fi  -g  g 

him 


